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In this work, Fermat’s principle is applied to derive a simple exact formula for refraction (reflection) in terms of the 
lengths of the incident and refracted rays. This formula is a nontrigonometric alternative to Snell’s law and is 
general for all optical surfaces. It is used to derive the paraxial optics equations in a more simple and direct 
way than that often used in the literature. It’s also applied to derive a new single, exact ray tracing formula 
for the nonparaxial refraction (reflection) at a single optical surface. The obtained formulas are used to develop 
a simple ray tracing procedure for meridional refraction through systems of spherical surfaces without the need to 
use any form of Snell’s law. Numerical examples are provided and discussed. © 2012 Optical Society of America 
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1. INTRODUCTION 


Fermat’s principle can be considered as the basis for geome- 
trical optics. Its original formulation states that for a single 
reflecting or refracting surface, the actual path that a light 
ray follows from one point to another via the surface is 
one for which the time required is a minimum. This particular 
formulation is called the principle of least time [1-3]. In the 
literature, this formulation is often used to derive Snell’s law, 
which is considered the main step in the derivations of para- 
xial and nonparaxial ray tracing formulas. In the present 
work, we apply Fermat’s principle in a different way to derive 
an exact, nontrigonometric formula that can, successfully re- 
place Snell’s law in many applications. This formula is written 
in terms of the lengths of the incident and refracted (reflected) 
rays and can be applied for all optical surfaces (spherical, as- 
pheric, planar, etc.). Using this formula, we can easily derive 
exact ray tracing formulas in both paraxial and nonparaxial 
regions. Section 4 is devoted to the application of this formula 
in deriving the well known paraxial (Gaussian) equations. The 
technique of this derivation was found to be simpler and more 
direct than the usual trigonometric technique based on Snell’s 
law. This may indicate the usefulness of this formula in teach- 
ing paraxial optics in the undergraduate level. 

In the literature, the exact ray tracing of nonparaxial refrac- 
tion (reflection) at a single optical surface is mainly divided 
into two classes: meridional and skew. In both cases the tra- 
cing of a single ray can be achieved by applying a procedure 
made up of a series of equations in which Snell’s law in its 
scalar or vector form is the cornerstone. There are several ex- 
cellent ray tracing procedures that can be used [4—7]. To the 
best of the author’s knowledge, there are no published works 
in which we can find a single, nontrigonometric formula for 
meridional refraction and reflection at a single optical surface. 
Such a single formula is missed in geometrical optics. An im- 
portant goal of this work is to derive this single, exact formula 
for cases of refraction [see Eq. (44)] and reflection [see 
Kgs. (64) and (65)]. Indeed, the reflection formula can be 
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directly derived from the refraction one by just setting 
n = -1. This single formula is expected to have some theore- 
tical and practical applications in optical design. Optics stu- 
dents and optical designers can easily use it to compute 
the position of the image formed by a paraxial or nonparaxial 
ray in a short time using a pocket calculator. Therefore, it can 
be considered an exact generalization for paraxial optics 
equations. Although the present paper is devoted to the case 
of spherical surfaces, the obtained exact formulas can also be 
applied for all other kinds of optical surfaces provided that 
both the radius and the center of curvature corresponding 
to the point of incidence are known. Another goal of this work 
is to derive a simple, exact ray tracing procedure for the case 
of a centered system of spherical surfaces (Section 5.C). This 
procedure consists of only five equations and doesn’t depend 
on any forms of Snell’s law, which makes it new and different 
from other known procedures. All the results obtained in this 
work are followed by numerical examples. 


2. APPLYING FERMAT’S PRINCIPLE TO A 
GENERAL REFRACTING (REFLECTING) 
SURFACE 


Figure 1 shows a general surface (spherical, aspheric, etc.) 
separating two media of refractive indices n, and n,. R is 
the radius of curvature corresponding to the point of inci- 
dence P. Let O and I be a point object and its image formed 
by the surface, respectively. The optical path lens of the ray 
from O to I is 


OPL = n,OP + PI. (1) 


Using the cosines law, we have 


A=OP= [R? + K? - 2RK cos w, (2) 


A'=PI = (R + K” + 2RK' cos wp. (3) 
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Fig. 1. Refraction by a single spherical surface. 


According to Fermat’s principle of least time, we have 


a _ 1 RK siny RK’ sin y _ 


PL 0. 4 
ae A 71 (4) 
Hence, we obtain 
A nk 
A aK @) 


In Fig. 1, we assumed that the refracting surface is converging 
(m1 >). A and A’ are the lengths of the incident and re- 
fracted rays, and K and K’ are the distances between the cen- 
ter of curvature C and both the object O and the image I/, 
respectively. A similar procedure can be followed for the case 
of a diverging surface, where we finally get 


A_ kK 


A “yk 


(6) 


In fact Eqs. (5) and (6) exactly describe the refraction of any 
meridional ray by an optical surface. The meaning of the ne- 
gative sign will be cleared in the sign convention (Section 3). 
For the case of a diverging mirror, the application of Fermat’s 
principle (following the same procedure) leads to the formula 


A K 
“=, 7 
A’ kK’ @ 
while, for a converging mirror, we get 
A K 
ae. 8 
A’ OK’ (8) 


It’s obvious that Eqs. (7) and (8) can be derived from Kas. (5) 
and (6), respectively, by setting m) = —n,. Therefore, we have 
an exact formula for reflection and refraction that can be writ- 
ten in the general form 


A NK 
= +7 _, 
A’ nk" O) 


It’s not necessary for K and K’ to be on the optical axis. It’s so 
easy to show that Eq. (9) remains correct if the line OJ meets 
the radius of curvature R (or its extension) at any arbitrary 
point C’. In this case, the image J will be formed at a point 
I located on the extension of OC’. This can be shown in Fig. 2, 
where K = OC’ and K' = IC’. In case of a planar refracting 


Fig. 2. Line OJ can meet the radius at any point C’. 


surface, the point C’ is any arbitrary point located on the nor- 
mal to the surface N, as in Fig. 3. 

For a centered system of spherical surfaces, and taking sign 
convention into account, we can apply Fermat’s principle for 
the mth surface (Fig. 4) to get the general formula 


Am = Nm-1Km (10) 
An NK m 


3. SIGN CONVENTION 


1. Light rays travel from left to right. 

2. The radius of curvature, R, is positive for a convergent 
(refracting/reflecting) surface and vice versa. 

3. The distance S, between the object and the vertex V of 
the surface is positive when the object is located to the left of 
V and vice versa. 

4. The distance S, between the image and the vertex V of 
the surface is positive if the image is located to the right (left) 
of the refracting (reflecting) surface and vice versa. 

5. The distance K between the object and the center of 
curvature C for a (reflecting/refracting) surface is positive if 
the object is to the left of C and vice versa. 

6. The distance K’ between the image and the center of 
curvature C for a (reflecting/refracting) surface is positive if 
the image is to the right of C and vice versa. 

7. The distance L between the object and the image is 
positive if the former is located to the left of the latter and 
vice versa. 

8. The length of a ray, A, is positive if the ray is real and 
vice versa (i.e., negative for an extension of a ray). 

9. Heights above the optical axis are positive and 
vice versa. 


Fig. 3. Refraction at a planar surface. The line OJ can meet the nor- 
mal N at any point C’. 
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Fig. 4. Refraction of a meridional ray through the last two surfaces of 
a centered system of m spherical surfaces. 


10. The angle that the ray makes with the optical axis is 
positive if the latter has to be rotated in the counterclockwise 
direction (through the acute angle) to coincide with the ray 
and vice versa. 


N.B. To check the validity of the above sign convention in 
various cases of reflection and refraction, we can use the fact 
that the separation between the object and the image is 
L=K+K’. ie, we can substitute for K’ in Eq. (10) 
by kK'=L-K. 

In the following sections, we apply Eq. (10) in the ray tra- 
cing of both paraxial and meridional refraction and reflection 
at spherical surfaces. 


4. PARAXIAL RAY TRACING 


In this section we shall see how the application of Eq. (10) in 
deriving the well known paraxial formulas is easier and more 
direct than the usual derivation procedure based on the 
trigonometric form of Snell’s law. 


A. Refraction at a Spherical Surface 

In Fig. 1, the refraction of the meridional ray OP at a spherical 
surface is exactly expressed by Eq. (5). In the paraxial region, 
the point of incidence P is very close to the vertex V. So, we 
can approximately write A = S,, A’ = S,, K = (S, +R) and 
K’ = (S, - R). Substituting into Eq. (5), we get 


N Ny Mm —-N 


Ss Ss, -F 


1) 


Equation (11) is the well known paraxial formula for a single 
spherical surface. If the refracting surface is divergent, a simi- 
lar procedure is applied, taking sign convention into account. 
In this case we find that A = S,, A' = -S,, K = S, +R, and 
k'=S,-R. (K' is a negative quantity, since J is located to 
the left of the center C.) Now, substituting into Eq. (6), we 
again get Eq. (11). 


B. Reflection at a Spherical Mirror 

Figure 5 shows the meridional reflection at a convex spherical 
mirror. Equation (5) can be used for the case of a convex 
spherical mirror by setting n,; = —ny. It is reduced to the 
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Fig. 5. Meridional reflection at a convex spherical mirror. 


paraxial form as P comes very close to V, where A= S,, 
A'=-S,,k =S,-R, K' = R-S, (according to the sign con- 
vention, K’ is negative, since I is to the left of C). Now, sub- 
stituting into Eq. (5), we get 


(12) 


Equation (12) is the expected paraxial form for a spherical 
mirror. For a concave mirror, a similar procedure is applied 
using Eq. (6) and we again get Eq. (12). 


C. Refraction by a Centered System of Spherical 
Surfaces 
Equation (10) exactly describes the refraction of a meridional 
ray by the mth surface of a centered system of spherical sur- 
faces. This equation can be reduced to a paraxial form by mak- 
ing use of Fig. 4. When P,, comes very close to the vertex 
Vm, Wwe can consider that A, = Sim_1-tn1, Am = Sm) 
Km = Sm_1—- tm-1 — Rm, Km = Sm — Rm. Then, by substituting 
into Eq. (10), we get 

Mm Nm-1 Nm — Mm-1 


+ = ; (13) 
Sm bn-1 = Sm-1 Rn 


This is a general formula for a system of m spherical surfaces 
in terms of S,, and S,,_;. It can be easily used to trace a ray 
refracted through a centered system of m spherical surfaces 
in the paraxial region. In this case, the distances S,, and S;,_) 
are measured to the vertices of the surfaces and not to the 
principal planes. It should be noted that for the first surface 
(m = 1), Eq. (13) will have the form of Eq. (11). In this case t, 
can be seen as the separation between the first surface 
(m = 1) and a preceding (imaginary) surface (m = 0) in con- 
tact to the first surface, ie., ¢, = 0. S, is the distance of the 
image formed by the imaginary (m = 0) surface. Because 
it’s located to the left of the surface, it’s considered a virtual 
image, and according to the sign convention S, must be a ne- 
gative quantity, and Eq. (11) is obtained. Using Eq. (10) and 
Fig. 4, the paraxial formula [Eq. (13)] can be rewritten in terms 
of K,, and K),, as follows In the paraxial region we have 


Am = Kn +n, (14) 


An = Ki, + Rn. (15) 


Now, substituting from Eqs. (14) and (15) into Eq. (10), we get 
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Nn he Nm-1 = Nm ~ %m-1 
Ul 
Kn Kn Rin 


(16) 


Both Kgs. (16) and (13) are equivalent. They are general for 
divergent and convergent surfaces. For the first surface 
(m = 1), the distance K, in Eq. (16) should be negative, since 
it represents the distance of a virtual image formed by an ima- 
ginary surface (m = 0) preceding the first surface. Since the 
image is located to the left of the curvature center, K, is ne- 
gative, and Eq. (16) can be written for the first surface as 


ny + Nn _ ny — No 
Ki Ki Ry, 


(17) 


Equation (17) is the paraxial form of a single spherical surface 
written in terms of K. 


D. Paraxial Thick Lens Equation 

Using Eq. (13), we can easily derive the classical paraxial for- 
mula for a thick lens as follows. For the second surface 
(m = 2) we have 


Me My — Ng-Ny 
Sp t-Sy R, 


(18) 


where Sz is the distance of the final image formed by the lens 
from the vertex Vo, and t, is the lens thickness (Fig. 6). Now, if 
the object is located at infinity, ie., Sy = ©, n» = Nz = 1 (air) 
and n; = n, in this case, the distance of the final image F' from 
the second surface will be the second focal length f,, of the 
lens measured from F' to V2. Thus, Eq. (18) can be written as 


1 n l-n 
+ = : 19 
Jy 4-81 Ry a 

1 n l-n 
= , 20 
Ff, Slt /S) Be aie 

where 
nR 
5, =, (21) 
- _ (n-)Dt, 

F=(1-t/S] E m= (22) 


Fig. 6. Refraction of a paraxial ray through a thick lens. 
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Thus Eq. (20) can be written as 


Eon 1l-n 
Se (23) 
Sy S, Re 


Substituting from Eq. (22) into Eq. (23), we get 


1 1 1 a 
Sp - Df, zl ue »* ae 


where f,, is the focal length measured from F to the second 
principal plane, while f,, is the focal length measured from F’ 
to the vertex Vo, fp =fy/E. 


E. Thick Lens Equation in Terms of K 

The concept of “principal plane” plays a very important role in 
paraxial optics. In Section 4.D, the appearance of the quantity 
E=([1-(n-1)t/nR,] was natural and we assumed that 
E/f, = 1/f,p. Now, using Eq. (16), we shall do the same steps 
applied in Section 4.D, and a new quantity C corresponding to 
E will appear. First, Eq. (16) is written for the second surface 
of a thick lens in the form 


Nn Mm _ NM-M 
— ; 25 
Ky + Ke Rp (25) 


In Fig. 7, we have Ky = K'{ - e, so we get 


Ng a _— ma -y 
-. i ’ 
Ki-e K, Ry 


(26) 


where e is the separation between the two curvature centers 
of the lens, 


e=t- [R, — Ry]. (27) 


The equation of the first surface of the lens is 
n/K, +1/K, = (n,- 1)/R, so if K, = ~, then 


1/K', = (m1 - Y/R. (28) 


Assuming that n, = Nz = 1 and n; = n, Eq. (26) becomes 


-l n l-n 
; —~+—= , (29 
Kill-e/Ki] kK; ( Ro ) ) 
Now let 


C =[l-e/K']. (30) 


Fig. 7. Refraction of a paraxial ray through a thick lens. 
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Then Eq. (29) becomes 


ate (Be (31) 


[Cay 


Substituting from Eqs. (28) and (82) into Eq. (81), we get 


C n-1 1 1 : e 
Ky =( n ) [p: zl aan »' Laer aa 


which is the paraxial thick lens formula written in terms of K%. 
Here, K, is the lens focal length measured from the focus F' to 
the center Cy of the second surface. Now, as was done in 
Section 4.D, C/K; can be replaced by 1/f, where f, is the 
focal length measured from the focus F to a point that may 
be called the “principal center” in analogy to the concept 
of the “principal plane.” So, Eq. (83) becomes 


1 n-1\[ 1 1 e 
Sc =( n Jz zl ge »* Se = 


The location of that assumed “principal center” was found to 
be exactly that of the principal plane, i.e., f, =f. This can be 
proved by substituting from Eqs. (27) and (29) into Eq. (84) 
where we get 1/f, = 1/f, which means that 


(ome ae 
Ky fo Sp 


(35) 


In Section 4, we formulated the paraxial formulas in terms of 
K, i.e., the distances measured to the center of curvature. This 
new formulation of the equations may be useful in some ap- 
plications in which the vertex position is hard to determine. It 
also facilitates the derivation of exact nonparaxial (meridio- 
nal) formulas in simple forms, as will be seen in the next 
section. 


5. MERIDIONAL RAY TRACING 
A. Refraction at a Single Spherical Surface 


1. Description 

In Fig. 1, the meridional ray OP is refracted at the spherical 
surface and the image I is formed at the optical axis. This 
refraction was exactly described by Eq. (5), which can be 
rewritten as 


A' = pk’, (36) 
nA 
Baye (37) 


where n = 7,/n,. Now, using the law of cosines we get 
A® = R? + K” + 2RK’ cos yw, (38) 


where 
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/p2 _ pd 
cos y = = (39) 


Substituting from Eqs. (86) and (39) into Eq. (38) and solving 
for K’, we get 


_ VER -W + VR WW 


K' . 
eo 


(40) 


Multiplying the nominator and denominator by the conjugate 
of the nominator, we get 


2 
K'= = = (41) 
of Pe BP RP 


and the final form of S; is 


R2 


a= Re 
, JPR —h2 — JR he 


(42) 


In Egs. (41) and (42) the positive sign is chosen if R is positive 
(.e., for a converging surface) and vice versa. Using Eq. (37) 
and Fig. 1, we can find / in the form 


2 
pak e+e VR =i] , (43) 
The sign of the small root VR? — h? is negative for a conver- 
ging surface and vice versa. Now, substituting for # from 
Eq. (43) into Eq. (42), we get 


S,=R+ ua 
nf A? +R? + 2K VR? —1?| - (2) - VRP =e 


(44) 


Equation (44) is a main goal of this work. It’s a single, exact. 
formula that can be used to calculate the position of the image 
formed by a nonparaxial (meridional) ray in both refraction 
and reflection. In the case of reflection, we just set n = -1, 
where Eq. (44) is reduced to a very simple form [see 
Kq. (65)], as we shall see in Section 5.B.1. 


2. Paraxial Approximation 

In the paraxial region, the Gaussian formula [Eq. (11)] can 
be directly obtained from Eq. (44) by setting h=0 and 
K=S,+R. Also, we can get the paraxial form of a single 
spherical surface written in terms of K’s as follows. As h tends 
to zero, Eq. (41) for a converging surface can be reduced to 
the form 


1 pl 
KR R ) 
and Eq. (43) reduces to 
nR 
=n-—. 46 
pan-Z (46) 


So, from Eqs. (45) and (46) we get 
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n 1 n-l 
Kin R 7) 


Equation (47) is the paraxial form for a single refracting sur- 
face. It’s general for both converging and diverging surfaces. 


3. Parallel Incidence 
If K = ~ (parallel incidence), then Kq. (44) directly reduces to 


R2 
+ ‘ 
Vn? R? -h? - SR? - he? 


S,;=R (48) 


In this special case, the longitudinal spherical aberration exact 
formula is 


R Rr 


us a1 Vm VE 


(49) 


4. Numerical Example 

In Fig. 1, let S,=10cm, R=8cm (ie, K= 
S,+R = 18 cm), h=3 cm, and n=2 (n=n,/n,. To find 
S, for the formed image J, we just substitute in Eq. (44) to 
easily get S; = 41.8499969 cm. We’ve chosen the positive sign 
for the second term of Kq. (44) and the negative sign for the 
root VR? — h?, since R is positive. For the case of a concave 
surface, we can consider the following numerical example: let 
S,=20cm, Rk = -8 cm(ie., K =S, +R = 12 cm),h = 3 cm, 
and n = 2. To find S; for the formed image, we just substitute 
in Eq. (44) to get S; = —-11.43958214 cm. We’ve chosen the 
negative sign for the second term of Eq. (44) and the positive 
sign for the root VR? — h?, since R is negative. 


5. Refraction at a Single Spherical Surface in Terms of 0 
In Fig. 1, we assume that the height h is unknown and the an- 
gle 0 between OP and the optical axis is already given. Now 
we have just two equations: Kqs. (42) and (43). To calculate S; 
from Eq. (42), we must first determine both h and /. To do so, 
the following steps are suggested. Assuming that 0, S,, R, and 
n (N = N2/n) are all given, we can write 


K=S,+R. (50) 


The normal 7, is given by 


r, =K sino, (51) 
A=K cos 0+ VR? -?r’, (52) 
h=Asin 60. (53) 


The sign of the root in Eq. (52) was found to be negative if R is 
positive (i.e., for a converging surface) and vice versa. Equa- 
tions (50)-(53) can be summarized into just two equations as 
follows. Substituting from Eq. (51) into Eq. (52), and then sub- 
stituting for A into Kq. (87), we get 


B=n cos 06+n,/ (R/K)? - sin? 6. (54) 
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Substituting for A into Eq. (53) and using Kq. (54), we get 


Kp 
n 


h= 


sin 0. (55) 


6. Summary 
In terms of 0, the meridional refraction at a spherical surface 
can be exactly traced by using the following three equations: 


B=n cos 0+n,/(R/K) - sin? 6, (56) 


h= AP sin 0, (57) 


+ RF 
VER —h2 - JR? 1h?” 


S,;=R (58) 


where n, So, R, 0, and K = S, + Rare all given. The sign of the 
root in (56) is opposite to the sign of the quantity (R/K). 


7. Numerical Example 

Using Fig. 1, let Sp = 9.25 cm, R= 19.5035816 cm (e., 
K =S) + R = 28.7535816 cm), 0 = 26.565°, and n=1.5. 
The value of K’ can be calculated by substituting in the three 
equations [Eqs. (56)-(58)]. The results are # = 0.5766518, 
h = 4,9434327 cm, and S; = —23.8971421 cm. 


B. Meridional Reflection at a Spherical Mirror 


1. Description 

In Fig. 8, the meridional ray OP is reflected at a concave mir- 
ror and a real image is formed at I. This reflection was exactly 
described by Eq. (8), which can be rewritten as 


A’ = BK, (59) 


where 


Al 2 
p=z=z w+ [K £ VR - 12). (60) 


Now, using the cosines law in Fig. 8, and following similar 


steps to those of Section 5.A.1, we finally get 


R2 
K=+ : ; (61) 
V PR -h? + JR? -h? 


Fig. 8. Reflection by a single spherical surface. 
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In Eq. (60), the sign of the root VR? — h? is positive for the 
convergent mirror and vice versa. In Eq. (61), the sign is po- 
sitive for a convergent mirror and vice versa. Equations (60) 
and (61) can be combined into a single equation for a sphe- 
rical mirror in the form 


K'=4 = 
fe A? +R? £ 2K VRP 1] n+ VP-P 


(62) 


Equation (62) can be directly derived from Eq. (44). This can 
be proved if we notice that for a spherical refracting surface 
K'=S, -R. So, Eq. (44) can have the form 


PR 
K'=F 
nf [K? + Re + 2K VR? = 1?| - (2) - VPA? 


(63) 


Therefore, by setting n = —1 in Eq. (63), we get Eq. (62). The 
larger root in the denominator of Eq. (62) can be easily rear- 
ranged so that Eq. (62) can have the final simple form 


a. et 2 h\? 
wien (5) (64) 


Which can also be rewritten in terms of Sy and S; as 


1 i 2 h\2 
= it ; 65 
Ros; Bay (3) Om) 


To the best of the author’s knowledge, Eqs. (44), (64), and (65) 
are all new and were never published before. They are also 
general for all surfaces (planar, spherical, aspheric, etc.) pro- 
vided that K and R (corresponding to the point of incidence) 
are determined. The meaning of K in the case of a planar sur- 
face was shown in Fig. 4, where K = OC’ and C’ is any arbi- 
trary point located on the normal to the surface N. In this case, 
R is the distance PC’. 


2. Paraxial Approximation 
In the paraxial region, h tends to zero. So, Eqs. (64) and (65) 
reduce to 


al, i. 2 
KtK 7B (66) 

1 1 2 
(67) 


B-S,  R=5, "8 


Equations (66) and (67) are equivalent to the well known para- 
xial form [Eq. (12)]. They are general for convex and concave 
mirrors. 


3. Parallel Incidence 
When the object is at infinity (K = ~), and for a given finite 
value of h, Eq. (64) reduces to 
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2 
Ks = (68) 
The image distance S, from the vertex is 
R 
S,=RF Je (69) 


In this case, the exact formula for longitudinal spherical aber- 
ration is 


R R 
L + 5 (70) 
2° 2./R? — nh? 


The sign of the second term in Eq. (69) is negative when R is 
positive and vice versa, while in Eqs. (68) and (70) it’s the 
same as that of R. 


4. Numerical Example 

Let’s find the exact image distance S, for the convex mirror of 
Fig. 5 if R = -7 cm (negative as the surface is divergent), 
S, = 15 cm, and h = 4 cm. In this case we directly substitute 
in Eq. (65) to get S; = —1.709481753 cm (a virtual image). 


5. Reflection at a Single Spherical Surface in Terms of 0 
In Fig. 8, if 6 is given and h is unknown, the same steps as in 
Section 5.A.5 are applied for a spherical mirror, where we 
finally obtain 


h= (cos 0+ \/(R/K)? - sin? 0)K sin 0, (71) 


= a h\2 
ete Ri! (3) (72) 


6. Numerical Example 

Let’s find K’ for the image formed by a spherical concave mir- 
ror (R = +6 cm) for a point object placed at a distance S, = 
13.2 cm from the mirror vertex at the optical axis and 
@ = 22.5°. We simply substitute in Eqs. (71) and (72) to geth = 
4.585262251 cm and K’ = 2.825833372 cm. 


C. Exact Procedure for Meridional Refraction through a 
Centered System of Spherical Surfaces 

1. Description 

In the literature, there are several procedures for meridional 
ray tracing [4—7]. In this section, another new procedure based 
on the exact formulas [Eqs. (41) and (43)] is presented. The 
procedure consists of five equations, and we derive it follow- 
ing the same steps as in Section 5.A.5. Let’s consider the last 
two surfaces of a centered system of m spherical surfaces as 
in Fig. 9. For the mth surface, we have 


1 = Ky sin On-1, (73) 


Am = Ky, COS On) + / R3, - 77. (74) 
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Using Eq. (73), we can write 


Am = Ky, COS Om-1 + V dies = Kin sin” Ini. (75) 


According to the sign convention, A,, is a negative length, 
since I,,_1P is not a real ray. Thus, Eq. (75) should be 
written as 


An — -(Kn cos Om-1 + V R?, a K?, sin? én) , (76) 


Rm = Am SIN On-1- (77) 
Using Eq. (87) we have 


Bn = NmA m 
m ae ee 


= (78) 
Nm 1K m 


Substituting from Kgs. (76) into (78), we get 


nN A 
Bn = ae (cos al = VBn/ Kn)? - sin” Ons); (79) 
m- 


and using Eggs. (76), (77), and (79) we get 


hn = (“e")KnPn sin Om-1 (80) 
Nn 


where @,,_; is the angle between P,,J,,_; and the optical axis. 
We can determine 9,,_; in terms of 6,,_9 using the formula 


K, 
sin Om. = Ges | sin Oyy_9- (81) 
To derive Eq. (81), we use Fig. 9 as follows: 
T2 = K 


mm 


n—1 sin Om-1, (82) 
rT, = Ky_1 SM On_2, (83) 


YT, SIN Bg — Nm_2 


(84) 


rT, SNA, Mn-1 


Fig. 9. Refraction through the last two surfaces of a centered system 
of spherical surfaces. 
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Dividing Eq. (82) by Eq. (83) and equating the result to 
Kq. (84), we get Eq. (81). In Fig. 9, we also find that 


Kn = Ka —tn-1 + Ry-1 = Ry. (85) 
2. Summary 


For the mth surface of the system, we have a procedure of five 
equations in the following sequence: 


: Nm 2K m1 \ 2 
0 =(eEe On-2> 86 
SIN On_1 (me) SIN On_2 ( ) 
Kn = Kit — tm-1 + Ry-1 7 Ry, (87) 


n ie 
Bm =-— (cos On-1 aE VBn/Kn)? - sin? Oma); 


Nm-1 
(88) 
Rm = (“)Knbn SiN On-1, (89) 
Nm 
2 
ne Rm (90) 


VBR, - he, ~ VR}, = hi, 


This procedure is different from other ray tracing procedures 
[4-7], since it doesn’t include any form of Snell’s law. The dis- 
tances are measured to the centers of curvature not to the 
vertices. 


3. Remarks 
1. For the final image, the value of the distance S,, is 
given by 


Sn = Kin + Rn. (91) 


2. The sign of the root in Eq. (88) was found to be the 
same as that of the quantity (R,,/Ky»). 

3. The sign of Eq. (90) is positive for the convergent sur- 
face and vice versa. 

4. For the first surface (m = 1), and according to the sign 
convention, the distance K, is negative, since we assume the 
existence of an imaginary surface (m=0) preceding to 
the first surface. In this case, the object O is assumed to be 
the image formed by that (m = 0) surface, and since this 
image is located to the left of the surface center of curvature, 
the distance K, is negative. 


4. Numerical Example 

In Fig. 10, we have a system of five spherical surfaces that 
form a final image J; for a real object O located on the optical 
axis at a distance S, from the vertex V, of the first surface. 
Let’s calculate the exact distance S; of the final image from 
the vertex V; of the fifth surface if we know that the incident 
ray inclination angle is 6, = 17.30972443, S, = 12 cm, R, = 
10 cm, Ry, = -8 cm, Rz = 12 cm, Ry = -10 cm, Rs = -8 cm, 
t, =5 ecm, & =5cm, fg =8 cm, =1em.n, = 1, n, = 1.2, 
ng = 1, nz = 1.5, ny = 1, N5 = 1.2. For each of the system 
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Fig. 10. Centered system of five spherical surfaces. 


surfaces (starting from the first one), we substitute into the 
five Eqs. [(86)-(90)] as follows. 


First surface: 
0, = 17.30972448° (given), K, = -22 cm (given), 
By = -0.733293273, 
h, = 4 cm, and K{, = -33.12062834. 
Second surface: 
6, = 9.479589747°, Ky = —20.12062834 cm, 
Py = -1.123526355, hy = 4.467759332 cm, 
Ks, = -55.03140479 cm. 
Third surface: 
@2 = 4.148784933°, K3 = —80.03140479 cm, 
P3 = -1.299007738, hz = 5.008144154 cm, 
K, = 37.33759883 cm. 
Fourth surface: 
@3 = —5.926743117°, K4 = 51.33759883 cm, 
Ba = —0.552990504, hy = 4.397067985 cm, 
Ki’, = 17.76842329 cm. 
Fifth surface: 
6; = —26.58363535°, kK; = 14.76842329 cm, 
Bs = —0.70684292, hs = 3.892892507 cm, 
Ki, = 22.16416375 cm. 
The distance S; of the final image J; is given 
by S; = K, +R; = 14.16416375 cm. 


Vol. 29, No. 12 / December 2012 / J. Opt. Soc. Am. A 2687 


6. CONCLUSION 


Fermat’s principle of least time was applied to derive an 
exact equation in terms of the lengths of the incident and 
refracted rays. This equation successfully replaced Snell’s 
law and was found to be useful in deriving the well known 
paraxial formulas in more simple steps. It was also applied 
to derive exact single formulas for refraction and reflection 
of nonparaxial (meridional) rays at spherical surfaces. The 
author thinks that these exact formulas are missed in geome- 
trical optics and that they may be useful in theoretical and 
practical applications. Based on the obtained equations, a sim- 
ple exact ray tracing procedure was derived that doesn’t de- 
pend on Snell’s law. In this work, the positions of the object 
and the images were determined relative to the curvature cen- 
ters, which facilitated the analysis and simplified the resulting 
formulas. 
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